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基于 特征 航 迹 简化 模型 的 中 心 航 迹 提取 研究 


(中 国民 航 大 学 天 津 市 空 管 运行 规划 与 安全 技术 重点 实验 室 ， 天津 300300) 


摘 要 : 针对 现 有 中 心 航 迹 提取 方法 存在 的 问题 : 进 近 转 谊 点 容易 缺失 ，2。 数 据 量 上 庞大、 运算 耗 时 ， 有 些 航 迹 点 没有 
代表 性 。 根 据 “以 直 代 曲 ” 思 想 简 化 航 迹 点 数据 ， 然 后 基于 特征 航 迹 的 航向 进行 聚 类 ， 求 得 平均 航向 并 得 到 中 心 航 迹 
定位 点 ， 最 后 使 用 B 样 条 曲线 拟 合 得 到 中 心 航 迹 。 仿 真 结果 表明 ， 使 用 新 方法 航 迹 点 数量 减 小 78.92%， 同 时 准确 地 得 
到 中 心 航 迹 ， 解 决 了 转弯 点 缺失 和 耗 时 的 问题 ， 验 证 了 新 方法 对 于 中 心 航 迹 的 提取 具有 准确 性 和 有 效 性 。 
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Research on central track extraction based on simplified track model 


Wang Lili, Peng Bo 
(Tianjin Key Laboratory of Atc Operation Planning & Safety Technology, Civil Aviation Administration of China, Tianjin 300300, 
China) 


Abstract: In view of the existing problems of the central track extraction method : 1. approach turn points are easy to be missing; 
2. the amount of data is huge and the computation time is time-consuming; 3. some track points are not representative. Firstly, 
simplify the track point data according to the idea of "replace curve by straight"in this paper, then cluster base on the 
characteristic track course, obtain the average heading and the center track locating point, finally, use the B spline curve to fit 
and obtain the central track, simulation results show that the number of track points decreases by 78.92% using the new method, 
at the same time, obtain the central track accurately, and solved the problem of turning point missing and time consuming, prove 
that the new method is accurate and effective for the extraction of the central track. 
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一 聚 类 算法 ， 认 为 其 在 航 迹 聚 类 中 存在 容易 陷入 局 部 最 优 搜索 的 
缺点 , 提出 了 采用 遗传 模拟 退火 算法 对 模糊 C 均值 聚 类 算法 进 
通过 分 析 航 空 器 的 历史 雷达 数据 ， 航 空 公 司 能 够 客观 地 评 。” 行 优化 ， 徐 涛 等 人 四 提出 了 基于 航 迹 间 面 积 加 权 的 航 迹 相似 怕 
价 飞 行 员 的 飞行 品质 、 飞 机 性 能 和 航班 的 运行 情况 ， 进 而 科学 。 ”度量 的 CURE 算法 ， 使 聚 类 效果 更 精确 ，Rehm 等 人 0 直接 利 
地 调整 飞行 计划 B-9; 对 管理 单位 来 说 , 可 以 利用 航 迹 数据 研究 。 用 已 有 航 迹 定义 其 相似 性 ， 针 对 同一 机 场 不 同 跑道 的 进 近 程序 
飞行 程序 的 设计 与 优化 以 及 噪声 预测 910。 然 而 空中 交通 流量 。 聚 类 。Gariel 等 人 0 进行 了 雷达 航 迹 聚 类 ， 并 对 航 迹 运行 进行 
与 日 俱 增 ， 航 空 器 的 航 迹 数据 也 呈 几 何 级 数 增长 ， 航 迹 数据 的 。 ”实时 监控 。Kenneth 等 人 63 针对 机 场 的 繁忙 区 域 ， 采 用 基于 网 
数量 十 分 庞大 ， 人 工 处 理 这 些 数据 显然 是 不 实际 ， 因 此 如 何 从 ， 格 的 聚 类 方法 对 不 同 机 场 的 航 迹 进行 分 析 ， 得 到 其 不 同 机 场 的 
海量 航 迹 数据 中 挖掘 出 有 价值 的 信息 越 来 越 受 到 民航 领域 研究 。” 进 离 场 中 心 航 迹 。 

者 的 重视 ， 国 内 外 很 多 学 者 对 航 迹 聚 类 进行 了 相关 研究 ， 赵 元 前 ， 在 对 航 迹 聚 类 的 研究 中 ， 进 近 转弯 点 中 心 航 迹 的 提 
棣 等 四 将 航 迹 视 为 以 机 场 雷达 为 原点 建立 的 立体 坐标 系 下 的 离 。，” 取 和 拟 合 是 研究 的 难点 。 在 对 进 近 航 迹 进行 聚 关 ， 提 取 中 心 航 
散 点 云集 合 ， 利 用 计算 几何 的 骨架 提取 算法 ， 求 得 代表 点 ， 经 。 迹 时 ， 最 终结 果 往 往 是 直线 航 段 能 够 得 到 中 心 航 迹 ， 而 转弯 航 
过 处 理 得 到 中 心 航 迹 并 识别 出 中 心 航 迹 的 特征 点 ， 王 洁 宁 等 人 段 则 容易 缺失 。 因 此 ， 如 何 对 进 近 转弯 点 的 历史 航 迹 进行 有 效 
外 分 析 了 国内 外 现 有 的 航 迹 聚 类 算法 ， 认 为 现 有 的 航 迹 聚 类 算 ”的 中 心 航 迹 提取 与 拟 合 具有 很 好 的 理论 价值 和 实际 意义 ， 本 文 
法 忽略 了 时 间 信 息 ， 而 航 迹 点 的 分 布 与 进 场 时 间 有 者 密切 的 联 。 ”提出 基于 特征 航 迹 简化 模型 的 中 心 航 迹 提取 方法 ， 根 据 “ 以 直 
系 ， 提 出 了 时 间 空间 的 聚 类 方法 ， 王 涛 波 等 人 9 研究 了 基本 的 。” 代 曲 ” 思 想 简化 航 迹 点 得 到 特征 点 ， 然 后 基于 特征 航 迹 的 航向 
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进行 聚 类 ， 对 聚 类 结果 中 的 每 一 类 别 的 特征 航 迹 求 得 平均 航向 
同时 得 到 中 心 航 迹 定位 点 ， 最 后 使 用 B 样 条 曲线 拟 合 得 到 平滑 zf 
的 中 心 航 迹 ， 本 文通 过 对 进 近 转弯 航 段 真 实 的 历史 航 迹 数据 进 | 从 多 
行 实验 分 析 ， 验 证 了 该 方法 的 准确 性 和 有 效 性 。 人 

200 上 上 
1 ”雷达 数据 处 理 和 进 近 过 程 描 述 -| 

本 文 要 处 理 的 数据 源 是 原始 二 次 雷达 数据 ， 首 先 对 原始 二 | 
次 雷达 数据 进行 解码 转换 成 全 文本 明码 类 型 ， 然 后 再 将 所 有 含 
有 有 用 信息 的 文本 条 目 按 关键 字 逐 一 捡 出 ， 并 将 所 需 字段 转换 io 航 段 三 
成 程序 需要 的 二 进 制 整数 或 浮 点 数 格式 ， 数 据 单位 保留 不 变 ， 而 向 余 全 全 20 疯 
最 后 整理 得 到 的 规范 航 迹 数据 形式 如 表 1 所 示 。 i 

图 1 进 近 航 段 
表 1 整理 得 到 的 二 次 雷达 数据 
时 间 。。 续 度 经度 > 航向 2 ”基于 特征 航 迹 的 简化 模型 
11:40:26 112.2890 32.1376 274 279 169.00 在 数学 和 物理 学 科 中 有 “以 直 代 曲 ”的 思想 ， 认 为 一 条 


线 可 以 用 在 曲线 上 取 点 连 成 的 折线 段 近似 表示 。 对 于 转弯 航 迹 
来 说 ， 也 可 以 利用 这 样 的 思想 来 做 处 理 ， 当 连续 的 某 几 段 航 迹 
走向 大 致 相同 时 ， 可 以 用 这 儿 段 连续 航 迹 的 起 点 和 终点 间 的 直 
虽然 可 以 从 该 表 的 雷达 数据 中 获取 航向 信息 ， 但 是 雷达 数 线段 来 近似 表示 。 本 文 把 这 里 的 起 点 和 终点 称 为 特征 点 ， 特 征 
据 每 4 秒 钟 采集 一 次 , 对 于 D 类 航空 器 , 转弯 速度 较 大 , 4s 时 ”点 间 的 直线 段 称 为 特征 航 迹 。 为 了 保证 准确 性 ， 使 特征 航 迹 尽 
间 过 长 ， 不 能 准确 地 刻画 出 航空 器 的 航向 变化 规律 ， 所 以 在 本 ”可 能 地 接近 原 航 迹 ， 本 文 希望 特征 点 是 表示 航 迹 走向 有 较 大 变 
文 第 2 节 中 重新 定义 了 特征 点 、 特 征 特 征 航 迹 以 及 航向 变化 等 ”化 的 点 。 假设 有 如 图 2 中 由 总 到 五 的 航 迹 点 构成 的 一 段 转弯 航 
几 个 概念 。 由 于 二 次 雷达 返回 的 是 经 纬度 信息 ， 只 有 转换 成 大 迹 , 可 以 看 出 人 和 五 航 迹 点 两 侧 航 迹 走向 明显 不 同 ， 且 差异 较 
地 华 标 才能 准确 的 反映 出 飞机 的 实际 运行 轨迹 ， 因 此 以 机 场 的 。 大 , 所 以 这 两 点 应 该 是 特征 点 。 由 CP 到 CP 组 成 的 特征 航 迹 可 
导航 台 为 坐标 原点 ， 建 立 直 角 坐 标 系 ， 将 返回 航 迹 的 经 纬度 数 。” ”以 体现 出 原 航 迹 的 趋势 。 

据 转 换 为 大 地 坐标 , 坐标 转换 得 到 直角 坐标 系 下 的 横 纵 坐标 ( 单 
位 为 干 米 ) 如 表 2 所 示 。 


11:40:30 112.2890 32.1376 274 226 183.00 


11:40:34 112.2870 32.1293 426 247 175.00 


表 2 经 纬度 坐标 转换 成 平面 坐标 
横 坐 标 纵 坐 标 


1 00173.93750 00244.25000 


2 00173.56250 00243.62500 


4 00173.25000 00243.00000 


图 2 以 直 代 曲 原理 图 


将 原始 数据 经 过 坐标 转换 处 理 后 用 MATLAB 作 图 得 到 该 
进 近 程序 在 实际 运行 中 的 情况 ， 如 图 1 所 示 ， 算 选 截取 其 中 三 在 所 有 航 迹 点 中 选择 特征 点 ， 而 特征 点 的 最 佳 选 择 要 遵循 
段 转弯 部 分 分 别 命名 为 航 段 1、 航 段 2 和 航 段 3。 该 程序 的 飞行 ” ”两 个 原则 :准确 性 原则 和 最 简 原则 。 准 确 性 原则 是 使 特征 航 迹 
过 程 大 致 如 下 : 脱离 航路 后 向 台 飞 行 ， 在 导航 台 上 方 转弯 ， 此 a a 而 最 简 原则 是 希望 特征 点 的 个 数 
为 航 段 1 的 部 分 ; 然后 直 飞 下 一 个 导航 台 (IAF)， 过 台 转 弯 到 这 两 个 原则 在 选取 特征 点 的 过 程 中 是 相互 对 立 的 ， 只 
下 一 个 导航 台 ， 然 后 再 过 台 转 弯 ， 此 为 航 段 2 的 部 分 ; 此 后 进 i i 
入 基线 转弯 程序 直至 落地 ， 此 为 航 段 三 的 部 分 。 由 得 到 的 数据 ”上 度 达 到 要 求 ， 为 了 兼顾 准确 性 原则 和 最 简 原则 ， 使 特征 点 的 选 
显示 , 平均 一 条 航 迹 一 个 转弯 过 程 大 约 有 20 个 点 迹 , 该 机 场 一 。 取 达 到 最 优 ， 本 文 引入 特征 航 迹 与 原 航 迹 的 差异 度量 值 一 一 特 
年 的 起 降 架 次 大 约 三 十 多 万 ， 也 就 是 说 一 年 的 进 近 转弯 数据 累 ”” 征 航 迹 差异 和 矩阵 。E 的 定义 为 原 航 迹 的 近似 长 度 与 特征 航 迹 
计 可 达 进 四 百 万 。 这 样 庞大 的 数据 量 用 计算 机 处 理 起 来 不 仅 耗 OO 
时 ， 而 且 需 要 占用 非常 大 的 内 存 ， 不 经 济 。 因 此 需要 一 种 方法 长度 ,特征 航 迹 表征 长 度 包含 三 个 部 分 , 分 别 为 : a) 特征 航 迹 
既 可 以 在 很 大 程度 上 减少 航 迹 数据 ,又 可 以 保证 航 迹 的 准确 性 。 Se 即 特 征 点 之 间 的 各 原 航 迹 点 到 特 
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征 航 迹 的 距离 之 和 ; c) 角度 差异 量 4d, ， 如 式 〈1) 所 示 。 
四 
1,90° < 0 <180° 
1 是 各 原 航 迹 段 长 度 ，09 是 各 原 航 迹 段 向 量 与 特征 航 迹 
向 量 间 的 夹 角 ， 由 式 〈1) 可 以 看 出 ， 当 0<90° 时 4d, 等 于 该 航 
段 两 端点 垂直 差异 量 之 差 ， 如 图 3 所 示 ， 原 航 迹 的 长 度 为 
d=|p,+p,|+|p,+p,l+|p,+p, 
特征 航 迹 表征 长 度 为 


(1) 


< 


其 中 : 
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以 往往 导致 最 后 的 中 心 定位 点 出 现 缺 失 或 偏离 轨迹 的 现象 ， 为 
解决 传统 聚 类 算法 得 到 的 中 心 定 位 点 不 准确 的 问题 ， 聚 类 后 每 
一 簇 都 代表 了 转弯 的 某 一 阶段 ， 对 每 一 簇 的 航 迹 求 平均 值 ， 再 
将 所 有 艇 的 数据 结合 起 来 可 以 反映 该 航 段 的 历史 航 迹 的 分 布 趋 
势 ， 也 就 是 中 心 航 迹 。 线 段 的 平均 值 可 以 用 某 一 条 直线 截取 线 
段 (下 称 截取 线 ), 计算 截 点 的 平均 值 。 然而 截取 线 不 能 随意 选 
取 ， 如 果 直 线 与 艇 内 航 迹 间 的 夹 角 偏 大 ， 那 么 航 迹 左右 两 边 则 
各 有 一 段 无 法 准确 计算 平均 值 ， 所 以 截取 线 应 该 与 大 部 分 航 迹 


=|p,+p,|+(di+d,)+|d,+(d, 一 d)+d,| 则 特征 航 迹 差 异 
量 为 已 =d -4d ,通过 控制 E 的 大 小 , 使 得 简化 选取 的 特征 点 结 
果 符 合 两 个 原则 。 综 上 ， 新 方法 在 解决 处 理 大 数据 耗 时 问题 的 
优势 就 体现 在 简化 模型 上 ， 该 简化 模型 可 以 在 初始 数据 处 理 时 
将 数据 精简 ， 从 而 减少 算法 耗 时 。 


Pa 
特征 航 迹 


图 3 特征 航 迹 差异 量 


3 ”基于 航向 的 聚 类 


上 文 提 到 传统 聚 类 算法 在 进 近 转弯 点 处 的 提取 效果 不 好 ， 
原因 如 下 : a) 按 照 基本 的 聚 类 方法 是 对 轨迹 分 布 在 数学 上 的 评 
价 指标 一般 考虑 距离 、 密 度 ) 进行 聚 类 ， 没 有 对 于 航空 器 的 
运行 过 程 进行 考虑 ， 最 终结 果 往 往 只 a ad 
体现 出 轨迹 的 线性 走向 ;b) 存 在 一 些 不 能 表征 航空 器 实际 运 
规律 的 航 迹 点 ， 要 合 去 此 类 航 迹 点 ， 因 此 本 文 提 er 
聚 类 算法 ， 通 过 该 聚 类 算法 可 以 将 转弯 过 程 体现 出 来 ， 航 空 器 
的 转弯 轨迹 在 空间 上 表现 的 特点 是 转弯 半径 大 , 通常 为 2~5 km， 
所 以 可 以 近似 认为 是 飞 几 个 直线 段 完 成 的 转弯 。 同 时 转弯 的 本 
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方向 垂直 。 因 此 ， 引 入 平均 向 量 表示 簇 内 大 部 分 航 迹 的 方向 。 
若 选取 的 截取 线 与 平均 向 量 相 垂 直 ， 则 能 够 准确 求 得 平均 航 迹 
的 长 度 可 以 得 到 增加 ， 中 心 航 迹 也 更 精确 。 在 处 理 的 时 候 需要 
将 航 迹 用 向 量 表示 ， 平 均 向 量 的 定义 由 式 (2) 给 出 。 


-V+ 十 .十 V+ 
a a 


2) 


n 

得 到 平均 向 量 后 ， 求 中 心 航 迹 只 需要 用 垂直 于 平均 向 量 的 
截取 线 扫 过 每 一 簇 中 的 航 段 ， 然 后 求 其 交点 的 平均 值 即 可 ， 因 
航 迹 点 间 的 线段 构成 。 在 求 平均 航 迹 时 ， 应 该 考虑 
每 一 簇 中 的 所 有 航 迹 点 ， 用 过 每 一 个 点 的 直线 与 航 迹 相交 的 点 
求 平均 值 ， 得 到 点 称 为 中 心 航 迹 定位 点 。 连 接 所 有 中 心 航 迹 定 
位 点 便 可 以 得 到 中 心 航 迹 。 但 并 不 是 所 有 航 迹 点 得 到 的 平均 值 
都 能 成 为 中 心 航 迹 定位 点 。 


A\ I、 | | 
(@ J 十 | | 


图 4 平均 航 迹 原理 图 


如 图 4 所 示 ， 假 设 某 一 篮 中 航 迹 如 图 分 布 ， 其 中 蓝 色 的 箭 
头 代 表 平 均 航向 ， 橘 黄色 线段 与 平均 航向 垂直 ， 红 色 轨 迹 表 示 
通过 计算 每 一 航 迹 点 的 平均 值得 到 的 平均 航 迹 。 可 以 看 出 ， 如 
a 
与 航 迹 交点 少 于 一 定数 量 时 ， 如 图 中 的 中 和 @@， 得 到 的 中 心 航 


质 是 航向 的 变化 ， 而 航向 的 变化 不 是 一 路 而 就 的 ， 而 是 逐渐 变 
化 的 过 程 。 考 虑 以 上 情况 ， 本 文 提 出 按照 转弯 航向 对 特征 航 迹 


进行 聚 类 ， 将 航向 差异 较 小 的 航 段 认为 是 同一 转弯 阶段 ， 通 过 
求 每 一 阶段 的 中 心 航 迹 从 而 最 后 还 原 出 整个 转弯 过 程 的 中 心 航 
迹 。 转 弯 阶 段 聚 类 的 具体 步骤 是 : 将 所 有 特征 航 迹 的 线段 转换 
为 向 量 形式 ,随机 选择 一 条 线段 ,计算 它 与 其 他 向 量 间 的 夹 角 ， 
通过 角度 阔 值 9 ， 得 到 与 其 方向 相似 的 向 量 作为 同一 簇 ， 从 而 
完成 聚 类 。 


4 中心 航 迹 定位 点 的 获取 


传统 的 聚 类 算法 将 所 有 的 航 迹 数据 进行 聚 类 ， 没 有 对 航 迹 
点 进行 分 类 进而 对 每 一 类 中 的 航 迹 点 聚 类 得 到 中 心 定位 点 ， 所 


Le 
点 阔 值 ， 只 有 当 交 点 个 数 大 于 某 个 值 时 ， 才 认为 此 时 得 到 的 中 
心 航 迹 定 位 点 是 可 取 的 。 甚 次 为 了 使 中 心 航 迹 更 平滑 ， 对 于 中 
心 航 迹 定 位 点 的 间隔 也 需要 进行 控制 。 如 图 中 的 和 @， 因 为 @) 
和 G@@ 相 隔 太 近 所 以 跳 过 不 取 ， 因 此 图 中 的 中 心 航 迹 应 该 为 由 定 
位 点 包 、 由 和 加 构成 。 


5 ” 样 条 曲线 拟 合 


由 上 节 程序 得 到 了 中 心 航 迹 定位 点 ， 将 中 心 航 迹 定位 点 连 
线 得 到 中 心 航 迹 。 为 了 使 最 后 得 到 的 中 心 航 迹 平滑 且 符 合 动力 


学 约束 ， 本 文选 择 使 用 三 次 B 样 条 曲线 进行 连接 。 样 条 曲线 在 
工程 设计 中 应 用 广泛 ， 是 经 过 一 系列 给 定点 的 光滑 曲线 ， 其 在 
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每 个 点 处 的 一 阶 和 二 阶 导 数 连续 ， 具 有 连续 、 曲 率 变 化 均匀 的 
特点 ， 因 此 符合 飞行 性 能 。 三 次 B 样 条 曲线 原理 如 下 : 
假设 在 区 间 (a，b) 上 给 定 一 个 分 割 : 
A:a=%<X<...<X < 大 = (3) 
在 〈a，b) 上 的 一 个 函数 SC9 称 为 差 值 三 次 样 条 函数 ， 如 
果 满 足下 列 条 件 : 
在 每 一 个 区 间 (Xt),i=12,.…n 内 S(X) 分 别 是 三 次 多 项 式 
函数 ， 在 节点 (i=42,.…,n 一 1) 处 成 立 : 
SX—0)=s" (x +0),k=0,1,2,...,n—1 (4) 
即 小 区 间 上 的 三 次 多 项 式 函 数 ， 在 拼接 点 处 过 具有 二 阶 连 
续 拼 接 。 


满足 差 值 条 件 y=S(X),i=0,1,...,n， 


王 莉 莉 ， 等 : 二 A 


其 他 参数 不 变 ， 取 选取 上 =0.20、0.21、 
验证 本 方法 的 优势 ， 使 
同 的 数据 进行 聚 0 


5s 左右 ; 为 了 
法 对 比 ， 如 


图 12 所 示 : 同样 选取 相 


法 耗 时 44s, 提取 得 到 17 个 中 心 航 迹 定 


> 人 
na 


天 
研究 


其 
迹 提 


0.22 时 ， 算 法 耗 时 都 在 
] KMEANS 算法 与 本 方 


次 


位 点 ,用 黑色 


性 
下 


从 图 中 可 以 看 出 KMEANS 算法 在 直 
是 在 转弯 点 处 出 现 了 缺失 现象 ， 尤 其 
离 中 心 航 迹 ， 得 不 到 理 

映 出 航空 器 实际 的 运行 规律 ， 
弯 点 、 解 决算 法 耗 时 以 及 获取 中 


说 明 传 
心 航 


选取 插值 间隔 为 0.5 时 得 到 33 个 控 


由 点 的 坐标 , 随后 将 j 


一 中 


是 第 10、 
想 的 中 心 航 迹 点 


随后 对 本 方法 得 到 的 17 个 中 心 航 迹 点 进行 B 样 条 拟 合 
制 点 , 表 4 中 列 出 了 部 分 控 
空 制 点 依次 连接 得 到 最 后 中 心 航 


线段 的 提取 效果 很 好 
11 个 点 ， 严 
旨 取 效果 ， 从 而 不 
统 聚 类 算法 在 提取 ; 
迹 的 问题 上 存在 不 足 。 


am mh > 


点 


| 


yc 
人 
芒 


< 


口 ， 


迹 如 图 13 


B 样 条 曲线 的 基 函 数 如 下 ; 
1 
G,,(1) = et +37 一 31 十 ])， 
G,(1) = lr 一 6 六 十 4)， 
te[0,1] G) 
G,,(?) = 本 + 3f +3t+1), 
G,,(1) = 2 
三 次 B 样 条 曲线 段 : 
1 4 1 0 p, 
_1 ， -3 0 3 0 
R&O= 训 可 3 -6 3 0| | (0) 
= “3 =3 小 | 及 
6 ”仿真 分 析 
本 文采 用 MATLAB 8.3 环境 编程 进行 仿真 ， 选 取 某 机 场 一 
套 进 近 程 序 的 转弯 段 航 段 三 一 周 的 航 迹 数据 ， 数 据 样 本 中 共有 
15341 个 航 迹 点 ， 算 法 中 需要 设置 的 三 个 参数 分 别 为 : a) 大 ， 
它 决定 了 特征 点 的 数量 和 特征 航 迹 近似 程度 ; b) sita ， 其 决定 


了 类 的 数量 和 聚 类 效果 ; c 


距离 和 交点 个 数 ， 决 定 了 拟 合 航 
的 天 值 ， 对 航 迹 点 进行 剔除 ， 使 得 航 迹 点 


原则 ， 得 到 K 值 、 
所 示 。 为 验证 
k=0.20、0.21、 


则 和 最 简 
结果 如 表 3 
数据 ， 选 取 

示 ， 其 中 蓝 色 空 心 圆 


涝 


过 于 集中 ， 


剔除 效果 ， 从 航 


代表 原始 航 迹 点 ， 攻 
则 下 的 航 迹 点 ， 从 图 中 可 以 看 出 , 选取 0.20 和 0.22 时 
存在 两 个 点 相 邻 的 情况 ， 而 选取 0.21 时 航 


) dist ， 它 控制 中 心 航 迹 定位 点 间 的 
迹 的 光滑 程度 。 首 先 选 取 不 同 
的 选取 符合 准确 性 原 
剔除 后 剩余 航 迹 点 和 剔除 率 的 不 同 
段 三 中 选取 一 条 航 迹 
0.22 时 的 剔除 效果 分 别 如 图 5~7 所 
色 实 心 点 代表 剔除 后 
的 航 迹 点 
迹 点 比较 


不 同 角度 阔 值 sita 进行 聚 关 
将 航 段 聚 成 八 类 ， 


聚 类 效果 很 好 ， 航 
的 间距 姜 值 dist ， 选 取 dist =0.5、 
所 示 , 当选 取 dist =1.0 时 航 迹 点 间距 


发 散 ,分布 较 均匀 ,符合 上 文 所 述 的 两 个 原则 , 所 以 选取 k=0.21 
进行 航 迹 点 剔除 ， 对 剔除 后 得 到 的 航 


迹 点 进行 聚 类 分 析 ， 选 取 
类 ， 如 图 8 所 示 当 选取 sita = 20" 时 ， 
段 无 缺失 ， 然 后 选取 不 同 
1.0 和 1.5 的 效果 图 如 图 9~11 
较 合适 , 航 迹 点 比较 分 散 ， 


所 示 , 图 中 蓝 色 空心 圆 代 表 17 个 中 心 航 迹 点 , 红色 十 字 代 表 使 
用 B 样 条 曲线 得 到 的 33 个 控制 点 。 
表 3 不 同 冰 值 k 得 到 的 结果 
闪 值 k 0.20 0.21 0.22 0.25 
剔除 后 航 迹 点 3327 3234 3167 2948 
剔除 率 78.31% 78.92% 79.36% 80.78% 
表 4 B 样 条 曲线 拟 合 得 到 的 部 分 控制 点 
坐标 点 序列 1 2 3 33 
控制 点 坐标 X 150.4979 150.3399 150.1492 147.6681 
控制 点 坐标 YY 139.1093 138.1118 137.0123 141.9310 
142 -一休 
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图 5 剔除 率 =0.20 的 效果 图 
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图 6 剔除 率 =0.21 的 效果 图 
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图 7 剔除 率 =0.22 的 效果 图 
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图 8 角度 取 20 时 聚 成 8 类 
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司 距 为 0.5 时 得 到 20 个 航 迹 中 心 点 


151 


图 10 
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间距 为 1.0 时 得 到 17 个 航 迹 中 心 点 
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图 11 间距 为 1.5 时 得 到 12 个 航 迹 中 心 点 
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图 12 KMEANS 算法 得 到 的 17 个 航 迹 中 心 点 
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图 13 最 后 拟 合 得 到 的 中 心 航 迹 


7 ”结束 语 


本 文 根 据 “ 以 直 代 曲 ” 思 想 提出 基于 特征 航 迹 简 化 模型 的 
中 心 航 迹 提取 方法 ， 分 别 经 历 了 原始 二 次 雷达 数据 转 码 、 原 始 
航 迹 点 经 纬度 坐标 转换 、 非 特征 航 迹 点 的 剔除 和 航 迹 点 简化 得 
到 特征 航 迹 点 ， 进 而 对 特征 航 迹 聚 类 分 析 ， 提 取得 到 中 心 航 迹 
点 定位 点 , 使 用 三 次 B 样 条 曲线 拟 合 得 到 平滑 的 符合 飞机 飞行 
性 能 的 中 心 航 迹 ， 通 过 对 真实 的 二 次 雷达 进 近 数据 仿真 同时 与 
其 他 算法 对 比分 析 ， 结 果 表 明 本 方法 在 转弯 点 处 得 到 的 中 心 航 
迹 点 无 缺失 现象 ， 在 处 理 相同 数据 时 新 方法 耗 时 更 少 ， 解 决 了 
现 有 方法 存在 的 问题 ， 验 证 了 新 方法 的 准确 性 和 有 效 性 ， 对 航 
迹 聚 类 分 析 尤 其 是 在 转弯 航 段 的 中 心 航 迹 提 取 方 面具 有 重要 的 
意义 与 价值 ， 下 一 步 将 围绕 终端 区 内 多 条 进 场 航 班 的 航 迹 数据 
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的 中 心 航 迹 提取 问题 进行 研究 。 
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